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Abstract 

We consider the problem of low-complexity channel estimation in <ligital ultra-wideband receivers. 
We extend some of our recent, sampling results for certain classes of non-bandlimihed signals ami develop 
several methods that take advantage of transform techniques to estimate channel parameters from a low- 
dimensional subspace of a received signal, that is. by sampling the signal below the NycniLst rate. By 
lowering the sampling rate, we reduce computational requirements compared to current digital solutions, 
allow for slower A/D converters and potentially significantly reduce power consumption of digital receivers. 
Our approach is particularly suitable for indoor wireless sensor networks, where low rates and low power 
consumption are required. One application to high-resolution acquisition in ultra- wideband localizers is 
also presented. ' 

Keywords 

UWB, multipath channel, singular value decomposition, anniliilating filters, waterfilling, filter banks, 
UWB localizers, 

I. Introduction 

Ultra-wideband (UWB) technology has recently received much attention duo to benefits 
of exceptionally large fractional bandwidth [2] [23], such as very fine time resolution for 
accurate ranging, imaging and multipath fading mitigation. UWB signals are generated by 
, driving an antenna with very short electrical pulses, typically on the order of a nanosecond, 
thus spreading the signal energy from near DC to a few gigahertz. Dun to simplified 
analog front-end and duty-cycling nature of pulse-based communication, ultra-wideband 
technology promises much lower power consumption and higher integration than might 
be conventionally feasible. Although the possibility of using extremely short, pulses for 
certain applications (such as ranging or imaging) has been investigated for at least two 
decades, there still remains a lot to be done for this technology to become pervasive using 
newly defined FCC regulations [5]. Some of the important issues that remain include 
low-cost and low-power designs and novel system architectures that allow for coexistence 
with other narrowband systems. 

One of the crucial tasks of an UWB "system, which imposes serious restrictions on the . 
system performance, is synchronization. There is a vast literature that has appeared 
recently [2J [4] [7] [9] [11] [23], addressing both algorithmic and implementation issues of 



several synchronization techniques, with a clear trend to eliminate as much as possible the 
necessity for analog components and perform all processing digitally 1 [4] [11]. 

While high-performance schemes have already been proposed for analog systems [6] [7], 
their application to digital- oriented solutions is still not feasible due to prohibitively high 
computational requirements. Furthermore, implementation of such techniques in digital 
systems would require very fast and expensive A/D converters (operating in the gigahertz 
range) and therefore high power consumption. And finally, they arc mostly based on 
exhaustive search and are inherently time consuming. 

We propose several frequency domain methods for channel estimation and synchro- 
nization in UWB systems, which use low-rate uniform sampling and well-developed algo- 
rithmic solutions. We extend some of oiir recent sampling results for certain classes of 
non-bandlimited signals [12] [22] to the problem of channel estimation in ultra- wideband 
channels and develop algorithms for estimating \inknown time delays and propagation 
coefficients by sampling a received signal below the Nyquist rate. Our approach takes 
advantage of the low-rank property and the algebraic structure of the data matrix under 
noiso-frco conditions. It leads to reduced computational requirements and faster acqui- 
sition compared to proposed digital techniques, thus allowing for a practical hardware 
implementation and low power consumption. It is particularly suitable in applications 
such as precise position location or ranging, and can also bo applied to other wideband 
systems, primarily for timing synchronisation and localization purposes. We also analyze 
the performances of our developed schemes in the presence of noise, and discuss possible 
modifications of the algorithms using a filter bank approach. 

II. Subspace Channel Estimation 

A. Problem Statement 

A number of propagation studies for ultra-wideband signals have been done, which take 
into account temporal pi-opcrtics of a channel or characterize a spatio-temporal channel 
response [2] [18]. A typical model for an impulse response of a multipath fading channel 

'Digital implementation offers many advantages, Mich as lower power, cheaper technology, full integration, 
robustness etc. 
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is given by 

L 

fc(*)-2>*(t-fc) (1) 

where ^ denotes a signal delay along the Z-th path while ai is a complex propagation coef- 
ficient which includes a channel attenuation and a phase offset along the Z-th path. While 
this model does not adequately reflect specific bandwidth-dependent effects, it provides a 
good characterization of the propagation channel used for diversity reception schemes (i.e. 
RAKE receivers). Equation (1) can be interpreted as saying that a received signal y(t) 
is made up of a weighted sum of attenuated and delayed replicas of a transmitted signal 
s(t): i.e. 

= + (2) 

where r)(t) denotes receiver noise. Note that the received signal y(t) has only 2L degrees of 
freedom, time delays t\ and propagation coefficients a/. Although these parameters can bo 
estimated using the time domain model (2). an efficient, closed-form solution is possible 
if we consider the problem in the frequency domain. 

Let Y(lj) denote the Fourier transform' of the received signal 

L 

Y^u) = £ ai S(u;)e- j « tl -h N{u) (3) . - 

where S(lj) and N(u)) are Fourier transforms of s(t) and rj(t) respectively. Clearly, spectral 
components are givon by a sum of complex exponentials, where the unknown time delays 
appeal* as complex frequencies while propagation coefficients appear as unknown weights. 
Thus the problem of estimating the channel parameters can be considered as a special case 
of classic harmonic retrieval problems.-, well-studied in spectral estimation literature [10]- 
[14]. There is a particularly attractive class of subspace or SVD-based algorithms, called 
super-resolution methods, which can resolve closely spaced sinusoids from a short record 
of noise-corrupted data. In [14] a State Space method is proposed to estimate parameters 
of superimposed complex exponentials in noise, while a similar ESPRIT algorithm is pre- 
sented in [13]. In [10] several SVD-basod techniques for estimating generalized eigenvalues 
of matrix pencils are addressed, which can be directly applied to our case, such as Direct 
Matrix Pencil Algorithm, Pro-ESPRIT and its improved version TLS-ESPRIT. Another* 



class of algorithms is based on the optimal maximum likelihood (ML) estimator [1]. While 
ML methods generally require Z-dimensional search and are computationally more consum- 
ing than the SVD-baacd algorithms, in some specific cases they can bo simplified to allow 
for a closed form minimization procedure. In [21] two such methods - Iterative Quadratic 
Maximum Likelihood (IQML) method and MODE algorithm are discussed. 

In the following, we develop an algorithm based on the State Space approach [14], that 
provides an elegant and numerically robust tool to estimate the channel parameters from 
a signal subspace. In particular, we show that it is possible to solve for all the parameters 
from a set of samples obtained by sampling a lowpass version of the received signal. In the 
noiseless case, the minimum required width of the frequency band is determined by the 
number of degrees of freedom per unit or time (which is typically on the order of MHz), 
while in the noisy case it is mainly determined by a received signal-to-noise ratio. We will 
show that even when the signal power is well below thermal noise level, it is possible to 
obtain high-resolution estimates from a low-dimensional subspace. While in all derivations 
we will bo considering the lowpass version of the signal, it is worth noting that this is not a 
requirement for the success of our methods, that is, the developed algorithms would work 
with a bandpass approximation of the signal as well; 

B. Algorithm outline 

Suppose that the received signal y(t) is filtered with an ideal lowpass filter H b = 
rect(-/ b ,/ fe ) and sampled uniformly at a sub-Nyquist rate R H > 2/ 6 . Prom the set of 
samples, we can compute the DFT coefficients Y[m) that correspond to the lowpass ver- 
sion of y(t). Denote by S[m] the DFT coefficients of the transmitted UWB pidse mid let 
Y s [m) = Y[m]/S[m], assuming that in the considered frequency band this division is not 
ill-conditioned. Define next a P x Q data matrix J as 



Y«[l] 

Y,m 



YAQ-i) \ 
YAQ] 



(4) 



\n[P-i) Y,[P] ... y;[p + q-2]/ 



If we let P,Q>L and z x = 6"***, the data matrix J can be written as J = UAV T , with 



matrices £/, A, and V defined as 
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V = 



A = diag ( a 2 a 2 
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2l ^2 -3 



1 \ 



^2 



„Q-1 



(5) 



(6) 



(7) 



The state space method is basfed on two properties of the data matrix J. The first one 
is that in the theoretical case of noiseless data, J has rank I*. This will allow us to reduce 
the noise level by approximating a noisy data matrix with an optimal, rank L matrix. 
The second one is the Vandermonde structure of U and V, that is, they both satisfy the 
so-called shift-invariant subspace property, 



U=U-$ and V = V-® 



(8) 



where $ is a diagonal matrix having zfs along the main diagonal, while (•) and (-) denote 
the operations of omitting the first and the last row of (•) respectively. Our algorithm can 
be thus summarized as follows: 

1. Compute the spectral coefficients Y[m] from the set of samples 



y n =</i 6 (c- nT„),y{t) >, n »!,...,£] 



(0) 



where T„ = \/B.„ and Li > 2L. 
2. Define a P x Q matrix J as 

/ nio] y„[i\ 

Ys[l\ Y„{2) 



J = 



YJQ-1] \ 



(10) 



Vn[P-l] Y,[P) ... Y S [P + Q- lj/ 



whoro Y s [m] = Y[m]/S[m) and P and Q satisfy the following relations: P,Q > L and 

3. Compute the singular value decomposition of J 

3 = U R h 3 V? + U n KV? (11) 

The first term consists of L principal components, where the columns of U s and V s cor- 
respond to principal left and right singular vectors of J respectively. The second term 
consists of remaining nonprincipals due to noise. 

4. Estimate the signal poles z t = e~ jwotl as generalized eigenvalues of the matrix pencil 

V*-zVs. - (12) 

In practice, the generalized eigenvalues of (12) are estimated by computing the eigenvalues 
of a matrix Z, defined as 

Z = V+-V a (13) 
where (•)**" denotes the pseudoinverse of (-). 

5. Find the propagation coefficients a t from the Vandermonde system 

L 

Note that we have converted the nonlinear estimation problem into the simpler problem 
of estimating the parameters of a linear model.. Nonlinearity is postponed for the last 
step, where in all spectral estimation methods the desired information is extracted from 
the estimated signal poles. Yet, the state space approach provides a lot of flexibility in 
model parameterization with respect to finite precision errors and sensitivity of estimated 
parameters to model mismatch [15]. 

Note that in the above algorithm we considered the lowpass approximation of the signal 
to estimate all the relevant parameters. In practice, however, it may be desirable to sample 
a bandpass version of the signal, while the best performance of our method is achieved 
if we choose a frequency band where the signal-to-noisc ratio is highest. It is also worth 
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noting that it suffices to use only a portion of tho signal bandwidth 2 and yet obtain high- 
resolution estimates of all the relevant parameters, due to the fact that the problem we 
consider belongs to the class of non-linear estimation problems. In Section III we will 
discuss possible improvements of the developed framework by sampling several bands at 
the same time, which can be of interest, for eocample. in the case of very low signal-to- 
noise ratios. While that approach yields better numerical performances due to the fact 
that almost all the signal energy can be used, it also results in increased computational 
requirements and complexity of the receiver. 

The performance of the State Space algorithm in the low signal-to-noise ratio regime can 
be further improved in the following way. Since in the noiseless case V* and V 3 span the 
same column space, we can extract L principal components from and V a by computing 
the joint SVD 

The signal poles z t = e~ ju3 ° tl are then estimated as generalized eigenvalues of an L x L 
matrix pencil V Rl — zV&. While this approach typically leads to bettor estimation accuracy 
than the original State Space algorithm, it requires that the SVD of the two matrices. J 
and [K, V s ] is computed. 

C. Computational Complexity and Alternative Solutions 

A major computational requirement of the developed algorithm is dominated by the 
singular value decomposition step 3 M which results in the overall computational order of 
0(M 3 ), where M = max(P,Q). On the other hand, we are often interested in estimating 
the parameters of only few strongest paths, therefore, computing tlie full SVD of the matrix 
J is not necessary* Alternatively, we can use some simpler methods from linear algebra 
[3] to find principal singular vectors, which have lower computational requirements and 
converge, very fast to the desired solution. We first give an outline of the Power Method 
that can bo used to compute only one dominant right (or left) singular vector of J. This 
can be of interest for initial synchronization, specifically in applications such as ranging 
or positioning. Later, we present its extended version applicable to the general case of 

3 Thc minimum required bandwidth depends on the signal-to-noisc ratio and the nuiubor of parameters to be 
estimated. In mast canes encountered in practice, it doesn't havo to exceed one third of the signal bandwidth. 
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estimating > 1 principal singular vectors. 
Power Method 

Consider a matrix F = JJ H of size Px?, and suppose that F is diagonizable, that is. 
A~ l FA = c1iAg(Ai,... t A P ) with A = [v 1: ... : v P ] and |A X | > |A 2 | > ...|A P |. Given y<°\ 
the power method produces a sequence of vectors y<*> in the following way: 

y<*> = z<*>/||z< fc >|| 2 (16) 

The method converges if A x is dominant and if y (0) has a component in the direction of the 
corresponding dominant eigenvector v x . It is easily verified that v u ... % v P are the right 
singular vectors of J. therefore, once the principal singular vector v x has been estimated., 
the signal pole z A corresponding to the strongest signal component is given by z A = v^Hi. 
A possible drawback of this method is that its convergence rate depends on |A 2 /Ai|. a 
quantity which may be close to 1 and thus cause slow convergence. Improved versions 
of the algorithm which overcome this problem are discussed in [3]. Note that the power 
method involves only simple matrix multiplications and has the computational order of 
£?(P 2 ). 

Orthogonal Iteration 

A straightforward generalization of the power method can be used to compute higher- 
dimensional invariant snbspaces, that is, to find M d > 1 dominant singular vectors. The 
method is typically referred to as Orthogonal Iteration or Subspace Iteration and can be 
summarized as follows. 

Given aPx M A matrix W< 0) , the method generates a sequence of matrices W (k) through 
the iteration 

Z {k) = FW< fc - x) (17) 
W <&) R (&) = Z (*) (q_jj factorization) (18) 

The computational complexity of the method is on the order of 0{P 2 M d ), and clearly, 
when M d = 1 the algorithm is equivalent to the power method. In practice, F is first 



reduced to upper Hessenberg form (that, is, F is aero below the first subdiagonal) and the 
method is implemented in a simpler way that avoids explicit Q - R factorization in each 
iteration. A more detailed discussion on this topic can be found in [3]. 

D. More Realistic Channel Models 

We will now extend our analysis to the more complex case of a channel that t,akcs into 
account certain bandwidth-dependent properties. Namely, as a result of the very large 
bandwidth of UWB signals, components propagating along different propagation patlis 
undergo different frequency selective distortion and a more realistic channel model for 
UWB systems is of the form 

1=1 

where pi{t) are different puke shapes that correspond to different propagation patlis. In 
this case, the DFT coefficients of the received signal are given by 

L 

Y[m] = S[m] flHflie-«*' + N[m] (20) 
z=i 

where Pi[m] are now unknown coefficients. One possible way to jointly estimate the 
relevant channel parameters is to use an approximation of the DFT coefficients Pi[m] with 
polynomials of degree d < R — 1, i.e. 

R-l 

PiH = J2^m r (21) 

Equation (20) now becomes 

L R-l 

Y[m) = S[m] ]T J2^' mre ~ iinUJOtt + ( 22 ) 
By denoting c z , r = a#v and Y H [m\ = Y[m]/S[m], we ol>tain 

Y*W = 22 ci>rm r e- jmuoti + *M (23) 

Similarly to the approach discussed in Section II-B, we can form the data matrix J as in 
(10). However, in this case, the state space method cannot be applied since the matrices 



U and V (made up of left and right singular vectors of J) no longer satisfy the shift- 
invariance property. We will therefore present an alternative method, based on annihilating 
filters, which is commonly used in error-correction coding and spectral estimation [19] [22]. 
The main idea behind this approach is to find the so-called annihilating filter H{z) = 
^f =0 ff[fcj 2 - fc that satisfies 

(if*y s )[m] = 0, Vm€Z (24) 
A \iseful property of such a filter is that it has multiple roots at z\ = e~~ jwatl [22], that is, 

L RL 

H{z) = " e^ tl z' l ) B = Y, H[k}z~ h (25) 

1=1 k-0 

In the following, we give an outline of the algorithm, while a more detailed analysis of the 
annihilating filters is presented in [19] and [22]. 

Annihilating filter method 

1. Find the coefficients H[k] of the annihilating filter 

L RL 
1=1 k=Q 



-A: 



(26) 



wluch satisfies 



In matrix form, the system (27) is eqi 
/ : : 



RL 

i/[m]*r fl [m] = ^ff[fc]i;[m--A:]=:0 1 Vm G Z 

k=Q 

livalent to 



(27) 



mi] 



Y S [RL] Y„[RL-l) 



\ 



Y„[-(RL - I)] 
Y,[-{RL - 2)] 

Y„[0) 



( #[0] \ 
H[l) 

\H[RL]J 



(28) 



Since there are RL + 1 unknown filter coefficients, we need at least RL + 1 equations, 
therefore, the number of DFT coefficients we have to compute is at least 2RL + 1. By 
setting H(0) = 1, at critical sampling (28) becomes 
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/ K[0] 



ni-i] 

K[0] 



Y»[-{RL - 1)] 




(29) 



This systems of equations is usually referred to as the high-order Yule- Walker system [19]. 

2. Find the values of tj by finding the roots of H(z). Recall that H(z) which satisfies (27) 
has multiple roots at z x — e~ i<aati , that is, 

L 

H(z) - |J(i - e-****- 1 )* (30) 
i=i 

At this point it is worth noting that while this is true in the noiseless case, in the presence 
of noise it is more desirable to estimate the time delays from L roots of H(z) which arc 
closest to the unit circle. 

3. Solve for the coefficients c/ >r by solving the linear system of at least RL equations in 
(23), 

L R-l 
i-1 r=0 



While the above described method uses standard computational procedures, for low values 
of signal-to-noise ratio it should be modified since any least-square procedure that deter- 
mines H directly from (29) has poor numerical precision. In practice, this problem can bo 
dealt with by using oversampling and an SVD step, where the matrix Y s is decomposed 
as 

Y s - U s AsK H + U n A n Vf (32) 

with the first term corresponding to the best (in the Probenius-horm sense) rank L ap- 
proximation of the matrix Y s . The filter coefficients H are then computed as 



H = -V K K- l U? 



(33) 



\Ys[RL]J 
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So far we have considered only the low-dimensional subspace of the received signal, 
obtained by filtering the signal with a bandpass filter. A possible modification of the 
sampling scheme, obtained by sampling several bands and estimating the channel from a 
larger subspace is considered in the next section. 

Ill, Filter banks and Waterfjlung approach 

As already pointed out, our framework is a frequency domain approach that allows for 
high-resolution estimation of channel parameters by exploiting a low-dimension property 
of the signal subspace. While in the noiseless case it is possible to estimate the parameters 
from any subspace of appropriate dimension, in the presence of noise the best performance 
of our algorithm is expected when the anti-aliasing filter is chosen such that a signal-to- 
noise ratio at its output is maximized. An alternative approach would be to use a filter 
bank at the receiver where each subbaiid is sampled at a rate determined by the filter 
bandwidth, while the sampled signals are then combined before running the estimation 
algorithm. That is, the set of coefficients Y s [m] is computed for each subband separately, 
and then combined to form the matrix J in (10). An.obvkms advantage of this approach 
is that almost all the signal energy is used, while the power is saved by running slower 
A/D converters in parallel. However, this results in increased computational requirements 
of the estimation algorithm. 

In the case when the channel parameters are estimated from adjacent subbands, the 
developed algorithms remain essentially the same, since J is a Hankel matrix made up of 
consecutive DFT coefficients A more interesting case, and perhaps more important 

in practice, is when the parameters are estimated from bands that are not necessarily 
adjacent. For example, if noise level in certain bands is relatively high or if some bands 
are subject to strong interference (such as interference from narrowband systems), it is 
desirable to estimate the channel by sampling only those bands where SNIR (signal-to- 
noise-pliLs-interference ratio) is relatively high. We will show that our developed algorithms 
can be adapted rather simply to handle this case. 

Assume for simplicity that we want to estimate the parameters by sampling only two 
bands B t and 5 2 . Let Y[m], m € [Mi, N x ] U [M 2l N 2 ], be the DFT coefficients computed 
from the set of samples in each band, and let Y^m] = Y[m]/S[m] assuming again that the 



above <livision is well-conditioned. Next construct a block-Hankcl data matrix J as 
/ Y a [M x ] Y,[M X + 1\ ... Y^ + Q-l] \ 



J- 



n[M 2 ] Y M [M 2 + 1) 



YAMx + Pl + Q-2] 
K[M 2 + Q-1] 



(34) 



\n[Afa + Pa-l] n[M 2 + P 2 ] ... K S [M 2 + P 2 + Q - 2]/ 
In the noiseless case, the matrix J can be written as J = UAV T , where £/, A, and V are 
now defined as 



z l 



*2 



z Vi 



22 2 3 



*2 



z 3 



\ ^M 2 +P 2 -l JWa+/%-l 
\ z l z 2 ~3 



~A/ 2 



JI/2+P2-1 



(35) 



A = cizagr ( Oi 02 03 



(36) 



V = 



/ 1 

21 



1 

22 



1 

^3 



~3 



1 \ 



<3-i 



(37) 



Clearly, the matrix V has the same Vandermoncle structure as in 7), meaning that the 
shift-invariance property (8) holds in this case as well, that is, 



V a. V.$ 



(38) 



where $ is the diagonal matrix having z^s along the main diagonal. Therefore, the signal 
poles can be estimated from the right singular vectors of J, using the algorithm described 
in Section II-B. However, it is interesting to note that the similar approach can be used to 
estimate the signal poles from the left singular vectors, which is of interest, for example, 



in the case when the number of rows in J is larger than the number of columns. Namely, 
the key is to observe the following property of the matrix U 

Z7 = £/.$ (39) 

where (0 stands for the operation of omitting the rows 1 and Pi + 1 of (■) , and similarly, 
(.) denotes the operation of omitting the rows Pi and Pi + P 2 of (•)- That is, the same 
shift-invariant subspace property can be exploited in this case as well, while the only 
modification in the developed algorithm is that the matrices U and |£ are constructed by 
removing the first and the last row respectively in each block of {/. The same idea can be 
generalized to the case when we sample multiple frequency bands. 

When there is additive noise, we should first extract the principal components by com- 
puting the singular value decomposition of J 

* 3 = U a h 3 V s H + U n KV? (40) 

and then estimate the signal poles z x = e"* 1 ' 0 *' as eigenvalues of a matrix defined as 

Z = V+-V 9 (41) 

or alternatively 

Z = U^-W 3 (42) 
It is clear that if we are interested in estimating only the dominant left/right singular 
vector (or a few of them), we can also use the power method or the method of orthogonal 
iteration, while the signal poles are then found from (41), or alternatively (42). Also note 
that the annihilating filter method from Section II-D remains the same, since it doesn't 
rely on shift-invariance property. Namely, the only requirement is that the rows of the 
matrix Y s in (29) are made up of consecutive coefficients Yi[ra], while Y s itself doesn't 
have to be a Hankel matrix. 

IV. Rapid Acquisition in UWB localizers 

One of the most interesting applications of om- framework can be found in ultra- 
wideband transceivers intended for low-rate, low-power indoor wireless systems, for ex- 
ample, in systems used for precise position location. Such UWB transceivers, called local- 
izers, have already been developed [6] and they use low duty-cycle episodic transmission 



of a coded sequence of impulses to ensure low-power operation and good performance 
in a multipath environment. Yet, the synchronization still presents a bottleneck in the 
transceiver design. The system developed by Aetherwvni [6] is analog and uses a com- 
plex method for sequence acquisition, based on the so-called Time Integrating Correlator 
(TIC), which implements a cascade of 32 correlators in analog circuitry and uses exhaus- 
tive search through all possible code positions. A drawback of this approach, apart from 
being analog, is that it is time consuming 3 since transmissions are spaced apart while the 
correlation window spans a small fraction of the sequence cycle time T ri/C j c . Furthermore, 
implementation of TIC alone can take up to 30% of the circuit area and tends to consume 
a major amount of the total power. An equivalent architectmc based on "mostly digital" 
conception is proposed in [4], which consists of a simple baseband analog section followed 
by a large digital processing back-end, where a cascade of 128 programmable matched 
filters is used for synchronization purposes. This approach, however, entails all the diffi- 
culties encountered in its analog counterpart. Besides, sampling Is achieved using an A/D 
converter designed to run at 2GHz, and developing alternative methods that would allow 
for simpler system architectiu-es is still an open problem. 

Our previous results can be directly extended to the case of timing synchronization 
(acquisition) in such systems by modeling the received signal y(t) as a convolution of L 
delayed, possibly different, impulses with a known coding sequence g(t), that is, 

2/(*) = 5>#i(t (43) 
The corresponding DFT coefficients Y[m] are given by 



where G[m] denote DFT coefficients of g(t). If we use the polynomial approximation (21) 
of the spectral coefficients P/[m], the total number of degrees of freedom per one cycle is 
2RL. Therefore, the signal parameters can be estimated using low-rate uniform sampling 
and the method we already described in Section II-D. This potentially leads to a few 
orders of magnitude faster acquisition, as will be shown later. 

a It takes Approximately 1000 cycles to acquire .synchronization. 
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(44) 



V. Performance Analysis 

Timing performance of our developed algorithms in UWB systems used for position 
location is illustrated in Figure 1 and Figure 3. All results are based on averages over 
500 trials, each witli a different realization of additive white Gaussian noise. Wc first 
considered the case where a coded sequence of 127 first-derivative Gaussian impulses is 
periodically transmitted, while a sequence duration spans approximately 20% of the cycle 
time T C|/c j c . Transmitted UWB pulse Is illustrated in Figure 1(a). while the transmitted 
sequence and a received signal within one cycle are shown in Figure 1(b). We considered 
a possibility of "imperceptible operation", that is, we assumed that the power spectral 
density of a received signal is below thermal noise level. Delay estimates for the case 
with eight propagation paths, but one dominant (with 70% of the total power), were 
obtained for various values of signal-to-noise ratio (E b /N 0 ) and the root-moan square 
errors (RMSE) of delay estimation for the dominant component is shown in Figure 1(c). 
We first considered the SVD-based approach and analyzed the performance of the method 
for different sampling rates. Our algorithm clearly yields highly accurate estimates (that 
is, with a sub-chip precision) for a wide range of SNRs. For example, with the sampling 
rate of one tenth the Nyquist rate (N A = iv" n /10) and E b /N 0 « -10d£, the time delay 
along the dominant path can be estimated with an RMSE of approximately one sample 4 , 
within only 30 cycles. Note that the main reason for using multiple cycles for estimation 
is to increase the effective processing gain by averaging samples taken over N c cycles. 
In the same figure, we show the timing estimation performance of the algorithm based 
on the power method and compare the results with the Cramer-Rao Bound (CRB) [20], 
which represents the lowest achievable RMS error for any unbiased estimator. The two 
presented methods yield essentially the same RMSE, and obviously, the performance of 
the algorithms improves as the sampling rate increases. Note that RAISE is very close to 
the CRB in the case of Nyquist-ratc sampling for all considered values of signal-to-noise 
ratio. 

We next analyze the case with two dominant signal components, each containing 35% 

4 Thc error is expressed in terms of number of samples taken at the Nyquist rate corresponding to the received 
UWB signal. 
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of the total power, and compare the performances of the SVD-based algorithm and its 
simplified version that uses Orthogonal Iteration. RMSE of time delay estimation along 
two dominant paths versus spacing between the two components (normalized to a pulse 
width w) is illustrated in Figure 1(d). While, in general, sampling at one tenth of the 
Nyquist rate is sufficient to obtain time delay estimates with an RMSE of less than 2 
samples, when the spacing between pulses is less than 2tu, higher sampling rates are 
required. For example, if the spacing between the two components is w, time delays can 
be estimated with an RMSE of approximately 1 sample, by sampling the received signal 
at one fifth the Nyquist rate. 

In Figure 2, we show the performance of our algorithm in the case of joint pulse shape 
and timing estimation. A coded sequence of 127 first-derivative Gaussian impulses is 
periodically transmitted over a channel with two propagation paths (again assuming that 
the sequence time spans 20% of the cycle time), while a received (single) UWB signal is 
made up of two pulses of different shapes, as illustrated in Fig. 2(a). The received noiseless 
and noisy UWB signals for SNR=0dB are shown in Fig. 2(b). The received signal was 
sampled uniformly over the entire cycle at one fifth the Nyqiust rate and averaged over 
N c = 60 cycles. We used the "annihilating filter method to estimate unknown time delays, 
by finding the roots of the annihilating filter, as well as different pulse shapes by polynomial 
fitting of the'DFT coefficients: As already pointed out in Section II-D, the signal poles 
(and thus the unknown time delays) can be estimated by choosing L = 2 zeros closest to 
the unit circle. Once the time delays of the pulses have been estimated, the corresponding 
pulse shapes are obtained by polynomial approximation of the DFT coefficients. In this 
case, we used a polynomial of order R = 20, which clearly yields a very good approximation 
of the received waveforms. 

In Figure 3, we show the performance of our method when the received pulse has no 
longer an ideal first-derivative Gaussian shape. We tested the method using a measm*ed 
waveform from the propagation experiment performed at the Berkeley Wireless Research 
Center, illustrated in Figure 3(a). We assumed that the received signal is sampled at 
one tenth (or alternatively one fifth) the Nyquist rate, while the frequency band used for 
estimation is shown in Figure 3(b). The RMSE of delay estimation for the dominant path 



vs. SNR is presented in Figure 3(c), and obviously, both the SVD approach and the power 
method lead to highly accurate estimates (again with a sub-chip precision). Note that 
with the approach based on matched filters [4], the resolution is Kmited by the bandwidth 
of the anti-aliasing filter, that is, 10 samples. Besides, the average acquisition time would 
be approximately five times longer, due to the fact that the method uses a combination 
of serial and parallel search. 

Another attractive feature of our scheme is good performance in the presence of a 
strong narrowband interference. Namely, UWB systems must contend with a variety of 
interfering signals from narrowband systems and building high-Q notch filters on chip 
would be technically very difficult. Since our method solves the estimation problem in the 
frequency domain, we can simply exclude spectral components that belong to frequency 
bands of interfering signals, while the state space algorithm can be easily modified to 
handle missing blocks of data [14]. The performance of our estimator in the prasence of a 
strong sinusoidal signal (SIR=-30dB) is illustrated in Figure 3(d). 

In all the methods we presented so far, estimation is carried out by sampling the received 
signal uniformly over T cyde . However, in the case when the sequence spans less than 10% 
of the cycle time, or if we are interested in estimated delays along several dominant paths, 
a better idea is to take a "multiresohition approach". That is, we eaai first obtain a rough 
estimate of the sequence timing, by taking uniform samples at a low rate over the entire 
cycle. Later, we perform precise delay estimation by increasing the sampling rate, but 
sampling the received signal only within a narrow window where the signal is present. 
This requires a "two-step" estimation, but at the same time can significantly reduce the 
computational requirements of the proposed schemes. 

VI. Conclusion 

We presented several methods for subspace parameter estimation in ldtra-wideband 
systems, which are based on our recent sampling results for certain classes of parametric 
non-bandlimited signals. Our approacli takes advantage of well-known spectral estimation 
techniques, requires much lower sampling rate and, therefore, lower complexity and power 
consumption compared to existing techniques. Besides, it leads to faster acquisition and 
allows for identification of more realistic channel models without resorting to complex 



algorithms. We specifically considered the application to indoor wireless sensor networks, 
where low rates and low power consumption are required. The developed algorithms can 
be also used in other UWB applications (mainly for synchronization) as well as in other 
wideband systems (such as wideband CDMA). 
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Claims 



1 . A method for estimating impulse response of a wideband communications 
channel as 

sum {k from 1 to L} a(k) delta(t-t(k)) 
from a received ultra-wideband signal, filtered with a lowpass/bandpass filter and 
sampled uniformly at a sub-Nyquist rate, comprising: 

(a) determining discrete-Fourier-transform coefficients y(j) and s(j) from the sampled 
signal and a transmitted ultra-wide-band pulse, respectively, 

(b) determining dominant singular vectors of a matrix having y(j-H-l)/s(j+i-l) as its 
ij-elements; 

(c) estimating signal poles from the dominant singular vectors and determining the 
time delays t(k) from the estimated poles; 

(d) determining the propagation coefficients a(k) from a Vandermonde system 
y(m)/s(m) = sum {k from 1 to L} a(k) y(k) A m (for m from 1 to L). 



2. The method of claim 1, wherein L is chosen as the nupiber of dominant singular 
vectors in step (b). 



3. The method of claim 1, wherein L is chosen as less than the number of dominant 
singular vectors in step (b). 




Pig. 1. Timing recovery in UWB systems (a) Transmitted first- derivative Gaussian impulse (b) Transmitted 
sequence of 127 UWB impulses (red) and received signal within one cycle (blue). Received signal-to-noise ratio 
is E b /No = -15dB. Coding is achieved uriUi a PN sequence of length 187, while tlie spacing between pulses 
(chips in the PN sequence) is assumed to be SO samples, (c) Root-mean square errors (RMSE) of delay 
estimation (in terms of number of samples) vs. SNR for the case with one dominant path and for different, 
values of the sampling rate. We assumed that the samples of the received signal are averaged over N c = 30 
cycles. Estimation performances of an SVD-based algorithm and its simplified version using the power method 
are compared with the Cramer-Rao Bound (d) RMSE of delay estimation for the case with two dominant 
signal components vs. spacing between received pulses. The spacing is normalized to the pulse width. 
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(d) 

Fig 2. Joint pulse shape and delay estimation (a) Received noiseless signal made up of tmo short pulses 
having different shapes (b) Received noisy signal (blue) and the noiseless pukes (blue), (c) Estimated shape of 
the. first pulse, (d) Estimated sliape of the second pulse: In both cases, we used a polynomial of order R - 20. 
to approximate the PFT coefficients of the received signal 




Fig. 3. (a) Received pulse waveform (including multipaths) and assumed transmitted pulse shape, (b) Normalized 
power spectral density (PSD) of the received pulse. We used one tenth or one fifth of the signal bandwidth 
for estimation, as indicated with dashed lines, (c) RMSB of delay estimation vs. SNR for the dominant 
path. We assumed that the samples of the received signal are averaged over N e = 30 cycles, and compared the 
performance of algorithms witfi SVD and tfie power method, (d) Estimation performance in the presence of a 
strong sinusoidal interference, SIR = — ZQdB. 
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